Simulation Optimization and Application of Shearer Strapdown Inertial Navigation System Modulation Scheme

The operating attitude of a shearer based on a three-dimensional (3D) space scale is the necessary basic information for realizing intelligent mining. Aiming to address the problem of the insufficient perception accuracy of shearers, in this paper, the rotation model of the actual turning mechanism of the strapdown inertial navigation system (SINS) of shearers is established, and the error propagation characteristics of different single-axis rotation modulation schemes are revealed. Through theory and simulation, the optimal rotation modulation scheme is determined to be the improved four-position turn–stop modulation with a rotation of <360°. The experiment shows that the 24 h positioning error of this scheme is 3.7 nmile, and the heading angle changes by 0.06°, which proves that this scheme can effectively improve the attitude perception accuracy of the inertial navigation system (INS). The field application of the shearer operating attitude perception based on this scheme shows that the positioning error after error compensation is 17% of that before compensation, and the heading angle error is 75% of that before compensation, which verifies that this scheme can significantly improve the accuracy of shearer operating attitude perception in field applications. This scheme can achieve higher precision perception accuracy based on SINS and has broad application prospects in the field of high-precision pose perception of coal mining machines, roadheaders, and other equipment.


Introduction
The operating posture of a shearer, including the heading angle, pitch angle, roll angle, and 3D acceleration of the body and cutting part, is based on the three degrees of freedom dynamic information on the 3D space scale, which can determine the position of the shearer and predict the 3D extension direction of the working face in the coal-bearing geological space. The operating posture is also the basic information for realizing the shearer's memory cutting, automatic height adjustment, intelligent control, intelligent perception, remote visual monitoring [1,2], etc. As a difficult problem of intelligent mining equipment, accurate perception of the shearer's running attitude can not only provide a way to detect and predict the production status of the intelligent chemical working face but also provide data feedback for the autonomous intelligent control of the shearer. Therefore, it is of great significance for the construction of an intelligent working face to realize the high-precision position and attitude state perception of the shearer based on rotary SINS.
Since the 1990s, traditional coal mining technology powers such as the United States, Germany, and Australia have begun to explore and research intelligent coal mining technology and equipment, and they initially realized remote sensing of coal mining equipment. In SINS mainly uses a carrier coordinate system (b system), navigation coordinate system (n system), rotation coordinate system (s system), and inertial coordinate system (i system) in the attitude calculation process. When the IMU remains stationary at the initial time, the s system, b system, and n system remain coincident, the origin is located in the center of the IMU, and the coordinate axis refers to the three orthogonal axes of the IMU.
IMU starts at angular speed ω and rotates around the celestial axis, and the shearer remains stationary during rotation, so at time t, the included angle between the x-axes of the two coordinate systems is β = ωt. At this time, the direction cosine array of the b system relative to the s system is as follows: Assuming that the ideal angular velocity and specific force of system b relative to system i are b ib ω and b f , respectively, the ideal output of IMU gyroscope and accelerometer is as follows: Under the excitation of angular velocity s is ω and specific force s f , the output of the gyroscope and accelerometer is as follows: s s a a n r ∇ is the accelerometer random drift error, g.
After expanding Equation (3), subtract the ideal output value of the inertial sensor and omit the second order small quantity. Then, the output error can be obtained as follows: SINS mainly uses a carrier coordinate system (b system), navigation coordinate system (n system), rotation coordinate system (s system), and inertial coordinate system (i system) in the attitude calculation process. When the IMU remains stationary at the initial time, the s system, b system, and n system remain coincident, the origin is located in the center of the IMU, and the coordinate axis refers to the three orthogonal axes of the IMU.
IMU starts at angular speed ω and rotates around the celestial axis, and the shearer remains stationary during rotation, so at time t, the included angle between the x-axes of the two coordinate systems is β = ωt. At this time, the direction cosine array of the b system relative to the s system is as follows: Assuming that the ideal angular velocity and specific force of system b relative to system i are ω b ib and f b , respectively, the ideal output of IMU gyroscope and accelerometer is as follows: where ω s bs = [0 0 ω] T and f bs = [0 0 0] T . Under the excitation of angular velocity ω s is and specific force f s , the output of the gyroscope and accelerometer is as follows: where ω s is is the gyro output value, deg/s; f s is the accelerometer output value, g; δK g is the gyro scale factor error, ppm; δK a is the accelerometer scale factor error, mA/g; δE g is the gyro installation error, deg; δE a is the accelerometer installation error, deg; ω s is is the ideal gyro input value, deg/s; f s is the accelerometer ideal input value, g; ε n is the gyro random drift error, deg; and ∇ n is the accelerometer random drift error, g.
After expanding Equation (3), subtract the ideal output value of the inertial sensor and omit the second order small quantity. Then, the output error can be obtained as follows: where C n s = C n b C b s and C b s = C s b T . According to Formula (4), after the IMU rotates 360 • relative to the shearer, the end value of the east and north gyro drift is 0; that is, the constant error of the inertial sensor perpendicular to the rotation axis is offset, which will not affect the accuracy of the navigation system.

Analysis of Rotary Modulation Scheme and Mechanism
3.1. Single-Axis Continuous Rotation Modulation 3.1.1. Single-Axis Continuous Rotation Modulation Scheme The s system, b system, and n system coincide at the initial time, and the single-axis continuous rotation modulation scheme is shown in Figure 2.

1.
The turning mechanism starts from point s at zero time, accelerates to point a through time t a with angular acceleration α, and the rotation angle during acceleration is θ a ; 2.
After the angular velocity reaches the set value ω, it will rotate continuously at a constant speed.
According to Formula (4), after the IMU rotates 360° relative to the shearer, the end value of the east and north gyro drift is 0; that is, the constant error of the inertial sensor perpendicular to the rotation axis is offset, which will not affect the accuracy of the navigation system.

Single-Axis Continuous Rotation Modulation Scheme
The s system, b system, and n system coincide at the initial time, and the single-axis continuous rotation modulation scheme is shown in Figure 2.
1. The turning mechanism starts from point s at zero time, accelerates to point a through time ta with angular acceleration α, and the rotation angle during acceleration is θa; 2. After the angular velocity reaches the set value ω, it will rotate continuously at a constant speed.

Modulation Mechanism Analysis of Single-Axis Continuous Rotation
The continuous rotation time of a single axis is 2 a n t t π ω = + , and the attitude angle error caused by gyro bias drift is given as follows:

Modulation Mechanism Analysis of Single-Axis Continuous Rotation
The continuous rotation time of a single axis is t = 2πn ω + t a , and the attitude angle error caused by gyro bias drift is given as follows: where ε s x is the gyro x-axis drift, deg/h; ε s y is the gyro y-axis drift, deg/h; and ε s z is the gyro z-axis drift, deg/h.
The attitude angle error caused by the gyro scale factor error is as follows: where K gx is the x-axis scale factor error of the gyroscope, ppm; K gy is the y-axis scale factor error of the gyroscope, ppm; K gz is the z-axis scale factor error of the gyroscope, ppm; ω ie is the earth angular velocity, deg/s; and ϕ is the latitude, deg. The attitude angle error caused by the gyro installation error is calculated as follows: Ψ x = ω ie cos ϕE gxy t a + ω ie sin ϕE gxz t a − ω ie sin ϕ αt 3 a 6 E gyz + E gxz sin 1 2 αt 2 a +E gyz cos 1 2 αt 2 a − 1 − ω ie cos ϕ sin 2θ a 4ω −E gxy + E gyx The single-axis continuous positive and negative rotation modulation scheme is shown in Figure 3.

1.
The turning mechanism starts from point s at zero time and rotates to point a with angular acceleration α to reach the set angular velocity ω. The rotation angle of the process is θ a and the duration is t a ; 2.
Then, rotate to point d at a constant speed of angular velocity ω. The rotation angle of this process is 2π − 2θ a and the duration is t c = (2π − 2θ a )/ω; 3.
Then, decelerate to the starting point s with the angular acceleration −α. The rotation angle is θ a and the residence time is t s ; 4.
Accelerate to point d with acceleration α again to reach the set angular velocity −ω; 5.
Then, rotate to point a at a constant speed with angular speed −ω; 6.
Then, rotate to the starting point s at a constant deceleration of angular acceleration α, and the dwell time is t s .

Single-Axis Continuous Positive and Negative Rotation Modulation Scheme
The single-axis continuous positive and negative rotation modulation scheme is shown in Figure 3.
1. The turning mechanism starts from point s at zero time and rotates to point a with angular acceleration α to reach the set angular velocity ω. The rotation angle of the process is θa and the duration is ta; 2. Then, rotate to point d at a constant speed of angular velocity ω. The rotation angle of this process is 2π − 2θa and the duration is tc = (2π − 2θa)/ω; 3. Then, decelerate to the starting point s with the angular acceleration −α. The rotation angle is θa and the residence time is ts; 4. Accelerate to point d with acceleration α again to reach the set angular velocity −ω; 5. Then, rotate to point a at a constant speed with angular speed −ω; 6. Then, rotate to the starting point s at a constant deceleration of angular acceleration α, and the dwell time is ts.

Modulation Mechanism Analysis of Single-Axis Continuous Positive and Negative Rotation
The modulation period of single-axis continuous positive and negative rotation is as follows: During this period, the attitude angle error caused by gyro bias drift is as follows:

Modulation Mechanism Analysis of Single-Axis Continuous Positive and Negative Rotation
The modulation period of single-axis continuous positive and negative rotation is as follows: During this period, the attitude angle error caused by gyro bias drift is as follows: The attitude angle error caused by the gyro scale factor error is as follows: where Ψ ry = 2K gy ω ie cos ϕ(t a + t s ) + (t c + t a )ω ie cos ϕ K gx + K gy + ω ie cos ϕ −K gx + K gy − sin 2θa ω + 2t a . The attitude angle error caused by the gyro installation error is as follows: The four-position turn-stop modulation scheme with a rotation of >360 • is shown in Figure 4. The residence positions and order of the turning mechanism are as follows: These positions are completely symmetrical with the acceleration and deceleration processes, and the residence time of each residence position is also the same.

1.
The turning mechanism starts from point A, with angular acceleration α accelerating the rotation to point a 2 to reach the set angular velocity ω. The process rotation angle is θ a ; 2.
Then, rotate to point b 1 at a constant speed with an angular velocity ω, and the rotation angle of the process is 0.5π − 2θ a ; 3.
Then, decelerate to point B with angular acceleration −α. The rotation angle is θ a and the residence time is t s ; 4.
Then, rotate to points C, D, and A in the same way, and the residence time at each point is still t s ; 5.
Then, reverse the rotation of points B, C, D, and A in the same way to form a complete modulation cycle.

Modulation Mechanism Analysis of Four-Position Turn-stop with Rotation >360°
The modulation period of this scheme is as follows: where ta is the acceleration time between two retention points, s; tc is the uniform rotation time between two retention points, s; and ts is the residence time at each retention point, s. The attitude angle error caused by gyro bias drift is as follows: The attitude angle error caused by the gyroscope scale factor error is as follows:

Modulation Mechanism Analysis of Four-Position Turn-Stop with Rotation >360 •
The modulation period of this scheme is as follows: where t a is the acceleration time between two retention points, s; t c is the uniform rotation time between two retention points, s; and t s is the residence time at each retention point, s. The attitude angle error caused by gyro bias drift is as follows: The attitude angle error caused by the gyroscope scale factor error is as follows: The attitude angle error caused by the gyroscope installation error is as follows: In summary, the four-position turn-stop modulation scheme with a rotation of >360 • has a very good error modulation effect, which can significantly improve the accuracy of SINS.

Four-Position Turn-stop Modulation Scheme with Rotation <360 •
The four-position turn-stop error modulation scheme with rotation <360 • is shown in Figure 5. The residence positions and order of the turning mechanism are as follows: The residence positions are completely symmetrical, but the rotation angle is not completely symmetrical.

1.
The turning mechanism starts from point A, with angular acceleration α accelerating the rotation to point a 1 to reach the set angular velocity ω. The process rotation angle is θ a ; 2.
Then, rotate to point c 1 at a constant speed, and the rotation angle is π − 2θ a and the rotation time is t c1 = (π − 2θ a )/ω; 3.
From point c 1 , decelerate uniformly to point C with angular acceleration −α, and the rotation angle during this process is θ a ; 4.
Accelerate the rotation from point C to point c 2 at the angular acceleration α to reach the set angular speed ω, and the rotation angle during this process is θ a ; 5.
Thereafter, it continuously rotates to point d 1 at a constant speed. The rotation angle in this process is 0.5π − 2θ a , and the rotation time is t c2 = (0.5π − 2θ a )/ω. It decelerates uniformly from point d 1 to point D at an angular acceleration of −α. The rotation angle in this process is θ a ; 6.
Then, reverse the rotation of points B and A in the same way to form a complete modulation cycle, and the residence time of points A, B, C, and D is t s . 5. Thereafter, it continuously rotates to point d1 at a constant speed. The rotation angle in this process is 0.5π − 2θa, and the rotation time is tc2 = (0.5π − 2θa)/ω. It decelerates uniformly from point d1 to point D at an angular acceleration of −α. The rotation angle in this process is θa; 6. Then, reverse the rotation of points B and A in the same way to form a complete modulation cycle, and the residence time of points A, B, C, and D is ts.

Modulation Mechanism Analysis of Four-Position Turn-stop with Rotation <360°
The modulation period of this scheme is as follows: The attitude angle error caused by gyro bias drift is as follows: The attitude angle error caused by the gyroscope scale factor error is as follows: ( ) The attitude angle error caused by the gyroscope installation error is as follows:   The modulation period of this scheme is as follows: The attitude angle error caused by gyro bias drift is as follows: The attitude angle error caused by the gyroscope scale factor error is as follows: The attitude angle error caused by the gyroscope installation error is as follows:

Improved Four-Position Turn-Stop Modulation with Rotation <360 •
According to the analysis in Section 3.4, it can be seen that the gyroscope ε is not completely compensated in the four-position turn-stop modulation scheme with a rotation of <360 • . From Equation (17), the ratio of residual zero-bias drift to total drift is as follows: Sensors 2023, 23, 4290 9 of 16 The above equation is further transformed into the following equation: According to Equation (21), there are two methods to reduce the residual drift ratio. The residual drift ratio can be reduced by extending the residence time at each position or changing the acceleration process of the turning mechanism. It is easier to extend the residence time at each position in practical applications, but it is more difficult to find the optimal acceleration process of the turning mechanism. The residual drift ratio is reduced by prolonging the residence time at each position.
2π cos(θ a ), c 2 = 2 t a t c1 + 3 4 ; then, Equation (21) can be simplified as follows: When the angular acceleration and angular velocity of the turning mechanism are constant, c 1 and c 2 are constants. According to Equation (22), prolonging the residence time at each position can significantly reduce the residual drift ratio. Suppose that the residence time of the turning mechanism in positions A and D is t s1 , and the residence time in positions B and C is t s2 ; then, in a modulation period, four times will the turning mechanism remain in four stay positions. The sum of the attitude error caused by the bias drift of the gyroscope is as follows: Taking the east direction as an example, the attitude angle error caused by the zero drift ε of the east gyroscope is as follows: cos(θ a ) + √ 2(t s1 − t s2 ) = 0, the following is true: Then, the east attitude angle error is zero, and the residual drift ratio is also zero, reaching the optimal compensation state.

Simulation Analysis
It can be seen from the theoretical analysis that the gyro bias drift is the main factor that causes a decline in the positioning accuracy. Therefore, the best modulation scheme can be selected by simulating and analyzing the modulation effect of the above scheme on the gyro bias drift. Matlab numerical simulation software was selected to simulate and analyze the above scheme. The parameter setting is shown in Table 1 and the simulation time was 72 h.
In each scheme, the positioning error caused by the constant drift error of the gyroscope in the direction of the vertical rotation axis is shown in Figure 6 can be selected by simulating and analyzing the modulation effect of the above scheme on the gyro bias drift. Matlab numerical simulation software was selected to simulate and analyze the above scheme. The parameter setting is shown in Table 1 and the simulation time was 72 h.
In each scheme, the positioning error caused by the constant drift error of the gyroscope in the direction of the vertical rotation axis is shown in Figure 6

Determination of Optimal Modulation Scheme
Based on the above theory and simulation analysis, it can be seen that the error of the single-axis continuous rotation modulation scheme is the minimum. However, under the

Determination of Optimal Modulation Scheme
Based on the above theory and simulation analysis, it can be seen that the error of the single-axis continuous rotation modulation scheme is the minimum. However, under the condition of a long endurance operation, the bias drift error will diverge over time, which will seriously affect the system accuracy, so this scheme should not be adopted.
The single-axis continuous positive and negative rotation scheme has the worst accuracy, which is inconsistent with the theoretical analysis. The reason is that the gyro bias drift error cannot be completely suppressed due to the asymmetry of the rotation process in the whole cycle, and the gyro bias drift error remains due to the alternating residence time of the positive and negative rotations, which ultimately affects the accuracy of the navigation system, so this scheme should not be used.
The four-position turn-stop scheme >360 • can completely modulate the gyro bias drift in the vertical axis direction within a complete modulation period, and this scheme will not produce other error-influencing factors that diverge over time. This scheme can be used as an alternative.
The inertial sensor drift error of the four-position turn-stop scheme with <360 • can be well compensated, but the gyro drift error cannot be effectively compensated, which will eventually affect the accuracy of the navigation system, so this scheme should not be used.
The improved four-position turn-stop scheme with <360 • cannot effectively compensate the gyro bias drift error, but after the theoretical adjustment of the four-position residence time of this scheme, the residual gyro bias drift error within the residence time just cancels out the residual gyro bias drift error in the process of uniform rotation, making the gyro bias drift error completely modulated in a complete modulation cycle. This scheme can also be used as an alternative.
In summary, the four-position turn-stop scheme with >360 • and the improved fourposition turn-stop scheme with <360 • are alternatives. The compensation effect of both schemes is good, and the system accuracy is basically equivalent. However, in practical applications, the four-position turn-stop scheme with >360 • requires the installation of conductive slip rings on the turning mechanism, which increases the cost of the navigation system. On the other hand, the conductive slip ring may introduce external errors when it works for a long voyage, thus affecting the system accuracy. The improved four-position turn-stop scheme with <360 • does not require additional conductive slip rings, so the system cost is relatively low, and it is relatively easy to achieve. Therefore, the improved four-position turn-stop scheme with <360 • is the best single-axis rotation modulation scheme in practical engineering applications.

Experimental Study on the Improved Four-Position Turn-Stop with <360 •
The experiment selected the SMT-1 three-axis turntable as the working platform for testing the inertial sensor, as shown in Figure 7. During the experiment, the platform provided the reference motion of the gyroscope and recorded and stored the actual angle change curve. The output of the gyroscope in the inertial measurement component was input to the signal acquisition computer through the RS-232 serial port, and data processing and comparison were carried out in the later stage. time of the positive and negative rotations, which ultimately affects the accuracy of the navigation system, so this scheme should not be used.
The four-position turn-stop scheme >360° can completely modulate the gyro bias drift in the vertical axis direction within a complete modulation period, and this scheme will not produce other error-influencing factors that diverge over time. This scheme can be used as an alternative.
The inertial sensor drift error of the four-position turn-stop scheme with <360° can be well compensated, but the gyro drift error cannot be effectively compensated, which will eventually affect the accuracy of the navigation system, so this scheme should not be used.
The improved four-position turn-stop scheme with <360° cannot effectively compensate the gyro bias drift error, but after the theoretical adjustment of the four-position residence time of this scheme, the residual gyro bias drift error within the residence time just cancels out the residual gyro bias drift error in the process of uniform rotation, making the gyro bias drift error completely modulated in a complete modulation cycle. This scheme can also be used as an alternative.
In summary, the four-position turn-stop scheme with >360° and the improved fourposition turn-stop scheme with <360° are alternatives. The compensation effect of both schemes is good, and the system accuracy is basically equivalent. However, in practical applications, the four-position turn-stop scheme with >360° requires the installation of conductive slip rings on the turning mechanism, which increases the cost of the navigation system. On the other hand, the conductive slip ring may introduce external errors when it works for a long voyage, thus affecting the system accuracy. The improved four-position turn-stop scheme with <360° does not require additional conductive slip rings, so the system cost is relatively low, and it is relatively easy to achieve. Therefore, the improved fourposition turn-stop scheme with <360° is the best single-axis rotation modulation scheme in practical engineering applications.

Experimental Study on the Improved Four-Position Turn-Stop with <360°
The experiment selected the SMT-1 three-axis turntable as the working platform for testing the inertial sensor, as shown in Figure 7. During the experiment, the platform provided the reference motion of the gyroscope and recorded and stored the actual angle change curve. The output of the gyroscope in the inertial measurement component was input to the signal acquisition computer through the RS-232 serial port, and data processing and comparison were carried out in the later stage. In order to verify the error modulation effect of the improved four-position turn-stop scheme with <360 • , the static performance index of the INS needed to be determined first. The error test results of the IMU 24 h static test under the condition that the turning mechanism does not rotate are shown in Figure 8. It can be seen from Figure 8a that when the IMU remained stationary, that is, when the system error was not modulated, the positioning error gradually accumulated within 12 h after the start of the experiment and slowly decreased between 12 and 17 h, but then the positioning error increased again, and the final positioning error was 18 nmile in 24 h. It can be seen from Figure 8b that the heading angle presented a similar sine wave change; it showed a trend of first increasing, then decreasing to zero, then increasing, then decreasing to zero, and then increasing within 24 h. The final heading angle changed to 0.02 • , and the amplitude of oscillation was 0.08 • .
IMU remained stationary, that is, when the system error was not modulated, the pos ing error gradually accumulated within 12 h after the start of the experiment and s decreased between 12 and 17 h, but then the positioning error increased again, an final positioning error was 18 nmile in 24 h. It can be seen from Figure 8b that the he angle presented a similar sine wave change; it showed a trend of first increasing decreasing to zero, then increasing, then decreasing to zero, and then increasing with h. The final heading angle changed to 0.02°, and the amplitude of oscillation was 0.0 (a) (b) The angular velocity and angular acceleration of the improved four-position stop scheme with <360° were selected as shown in Table 1. The residence time at each position was 100 s, and the total test time was 24 h. This experimental scheme can pletely modulate the influence of inertial sensor bias error on the final error of the na tion system. The final positioning error and heading angle change in this schem shown in Figure 9. It can be seen from Figure 9a that the positioning error gradual creased within 15 h after the start of the experiment, and then decreased gradually final 24 h positioning error was 3.7 nmile, which is more than four times higher tha static experimental positioning accuracy. It can be seen from Figure 9b that the chan heading angle within 24 h showed a trend of first increasing, then decreasing to zero increasing, and then decreasing. The final heading angle changed by 0.06°, which is less than the static error. In addition, due to the superposition of the indexing mecha and the inertial navigation itself, the experimental system had significant heading noise. The angular velocity and angular acceleration of the improved four-position turnstop scheme with <360 • were selected as shown in Table 1. The residence time at each stop position was 100 s, and the total test time was 24 h. This experimental scheme can completely modulate the influence of inertial sensor bias error on the final error of the navigation system. The final positioning error and heading angle change in this scheme are shown in Figure 9. It can be seen from Figure 9a that the positioning error gradually increased within 15 h after the start of the experiment, and then decreased gradually. The final 24 h positioning error was 3.7 nmile, which is more than four times higher than the static experimental positioning accuracy. It can be seen from Figure 9b that the change in heading angle within 24 h showed a trend of first increasing, then decreasing to zero, then increasing, and then decreasing. The final heading angle changed by 0.06 • , which is 0.02 • less than the static error. In addition, due to the superposition of the indexing mechanism and the inertial navigation itself, the experimental system had significant heading angle noise.

Analysis on the Compensation Effect of Shearer Operating Attitude Perception Error
The SINS of the coal shearer in the No.211 working face of a mine is mainly use sense the operating posture of the coal shearer in real time and guide the autom straightening of the scraper conveyor, as shown in Figure 10. It is regarded as a horizo coal seam, since the average dip angle of the coal seam is less than 1°. The shape o scraper conveyor depends on the output coordinates of the shearer SINS in X and Y a First, the compensation effect of the positioning error of the shearer SINS was investig by analyzing the shape of the scraper conveyor.

Analysis on the Compensation Effect of Shearer Operating Attitude Perception Error
The SINS of the coal shearer in the No.211 working face of a mine is mainly used to sense the operating posture of the coal shearer in real time and guide the automatic straightening of the scraper conveyor, as shown in Figure 10. It is regarded as a horizontal coal seam, since the average dip angle of the coal seam is less than 1 • . The shape of the scraper conveyor depends on the output coordinates of the shearer SINS in X and Y axes. First, the compensation effect of the positioning error of the shearer SINS was investigated by analyzing the shape of the scraper conveyor.
The SINS of the coal shearer in the No.211 working face of a mine is mainly used to sense the operating posture of the coal shearer in real time and guide the automatic straightening of the scraper conveyor, as shown in Figure 10. It is regarded as a horizontal coal seam, since the average dip angle of the coal seam is less than 1°. The shape of the scraper conveyor depends on the output coordinates of the shearer SINS in X and Y axes. First, the compensation effect of the positioning error of the shearer SINS was investigated by analyzing the shape of the scraper conveyor. The automatic straightening interface of the scraper conveyor is shown in Figure 10b. The yellow line in the figure is the reference line, the green line is the actual path, and the red line is the straightening target. The output positioning information of the SINS shows that most of the scraper conveyor deviated from the reference line by about 800 mm, which is consistent with the cutting depth of the shearer by 0.8 m. However, the scraper conveyor between No. 118-150 supports presented an obvious bending shape, with an average deviation of 130 mm from the cutting depth of the shearer. The deviation between the No. 122 support and the reference line was the largest, reaching 1102 mm, and the deviation from the cutting depth of the shearer was about 300 mm. Figure 11 shows the operation track of the shearer before and after error compensation. After the error compensation, the average deviation of the scraper conveyor between No.118-150 supports was 822 mm from the reference line, and the average deviation of the shearer section depth was 22 mm. The straightening error of the scraper conveyor after the error compensation, The automatic straightening interface of the scraper conveyor is shown in Figure 10b. The yellow line in the figure is the reference line, the green line is the actual path, and the red line is the straightening target. The output positioning information of the SINS shows that most of the scraper conveyor deviated from the reference line by about 800 mm, which is consistent with the cutting depth of the shearer by 0.8 m. However, the scraper conveyor between No. 118-150 supports presented an obvious bending shape, with an average deviation of 130 mm from the cutting depth of the shearer. The deviation between the No. 122 support and the reference line was the largest, reaching 1102 mm, and the deviation from the cutting depth of the shearer was about 300 mm. Figure 11 shows the operation track of the shearer before and after error compensation. After the error compensation, the average deviation of the scraper conveyor between No. 118-150 supports was 822 mm from the reference line, and the average deviation of the shearer section depth was 22 mm. The straightening error of the scraper conveyor after the error compensation, that is, the plane positioning error, was about 17% of that before the compensation. The positioning accuracy of the shearer SINS was significantly improved. that is, the plane positioning error, was about 17% of that before the compensation. The positioning accuracy of the shearer SINS was significantly improved. Since the roll angle and pitch angle of the shearer fluctuate relatively little, the heading angle (which has a larger fluctuation) was selected to study the attitude angle error compensation effect of the SINS. After removing invalid data, the output fluctuation of 3571 heading angles during the period 00:00:01-07:16:26 on 22 July 2020 is shown in Figure  12. Before error compensation, the heading angle fluctuated significantly in this period, with a maximum fluctuation value of 5.03° and a final cumulative error of 0.92°. After error compensation, the characteristics of heading angle fluctuation were significantly reduced: the maximum fluctuation value was 2.56° and the final cumulative error was 0.69°. Since the roll angle and pitch angle of the shearer fluctuate relatively little, the heading angle (which has a larger fluctuation) was selected to study the attitude angle error compensation effect of the SINS. After removing invalid data, the output fluctuation of 3571 heading angles during the period 00:00:01-07:16:26 on 22 July 2020 is shown in Figure 12. Before error compensation, the heading angle fluctuated significantly in this period, with a maximum fluctuation value of 5.03 • and a final cumulative error of 0.92 • . After error compensation, the characteristics of heading angle fluctuation were significantly reduced: the maximum fluctuation value was 2.56 • and the final cumulative error was 0.69 • . After compensation, the heading angle error was 75% of that before compensation, and the attitude angle perception accuracy of the shearer SINS was also significantly improved.
Since the roll angle and pitch angle of the shearer fluctuate relatively little, the heading angle (which has a larger fluctuation) was selected to study the attitude angle error compensation effect of the SINS. After removing invalid data, the output fluctuation of 3571 heading angles during the period 00:00:01-07:16:26 on 22 July 2020 is shown in Figure  12. Before error compensation, the heading angle fluctuated significantly in this period, with a maximum fluctuation value of 5.03° and a final cumulative error of 0.92°. After error compensation, the characteristics of heading angle fluctuation were significantly reduced: the maximum fluctuation value was 2.56° and the final cumulative error was 0.69°. After compensation, the heading angle error was 75% of that before compensation, and the attitude angle perception accuracy of the shearer SINS was also significantly improved.

Conclusions
(1) This paper theoretically analyzed the propagation rules of inertial sensor drift error, scale factor error, and installation error in single-axis continuous rotation, single-axis continuous positive and negative rotation, >360° four-position stop, <360° four-position stop, and the <360° improved four-position stop error modulation scheme. It was proved theoretically that the <360° improved four-position stop scheme can eliminate

Conclusions
(1) This paper theoretically analyzed the propagation rules of inertial sensor drift error, scale factor error, and installation error in single-axis continuous rotation, singleaxis continuous positive and negative rotation, >360 • four-position stop, <360 • fourposition stop, and the <360 • improved four-position stop error modulation scheme. It was proved theoretically that the <360 • improved four-position stop scheme can eliminate the influence of gyro drift in the vertical rotation axis direction. Moreover, it achieves a relatively good error self-compensation effect. (2) In this paper, five single-axis rotation error modulation schemes were studied through simulation analysis. The research shows that the positioning error of the improved four-position turn stop scheme <360 • was about 0.01 nmile within 72 h of simulation time. After a comprehensive analysis, this scheme was determined to be the best rotation modulation scheme in practical engineering applications. Based on the optimal rotation modulation scheme, a single-axis rotation error modulation experiment was carried out, which verified that the scheme can effectively improve the sensing accuracy of SINS. (3) The field application of shearer operation attitude perception was analyzed and studied. The field application showed that the plane positioning error after error compensation was about 17% of that before compensation, and the heading angle error was 75% of that before compensation, which verifies the effectiveness of the error compensation algorithm in this paper, and the shearer operation attitude perception accuracy has been significantly improved. (4) Although the single-axis rotation error modulation of the SINS of the shearer was studied in this paper, there is a lack of specific research on the rotation error effect of the turning mechanism in the vibration environment, which will affect the system error self-compensation effect of the single-axis rotation modulation IMU applied to the long flight and the selection of the optimal attitude solution structure of the shearer. Therefore, the rotation error effect in the vibration environment needs to be further explored.
Author Contributions: G.W. and Y.S. conceived and designed the experiments; G.W. and X.F. wrote the paper; X.F. and M.L. reviewed and revised the paper; N.C. and J.L. performed the experiments; and F.Q. analyzed the data. All authors have read and agreed to the published version of the manuscript.
Funding: This work was supported by the National Natural Science Foundation of China (nos. 52104167, 51874276, and 52004273) and the Fundamental Research Funds for the Central Universities (no. 2020ZDPY0209).
Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.
Data Availability Statement: All data and code used or analyzed in this study are available from the corresponding author on reasonable request.